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ABSTRACT. Airborne and spaceborne altimeters provide measurements of sea-ice elevation, from 
wbicb sea-ice freeboard and thickness may be derived. Observations of tbe Arctic ice pack by satellite 
altimeters indicate a significant decline in ice thickness, and volume, over the last decade. NASA's Ice, 
Cloud and land Elevation Satellite-2 (ICESat-2) is a next-generation laser altimeter designed to continue 
key sea-ice observations through the end of this decade. An airborne simulator for ICESat-2, the Multiple 
Altimeter Beam Experimental Lidar (MABEL), has been deployed to gather pre-launch data for mission 
development. We present an analysis of MABEL data gathered over sea ice in the Greenland Sea and 
assess the capabilities of photon-counting techniques for sea-ice freeboard retrieval. We compare 
freeboard estimates in the marginal ice zone derived from MABEL photon-counting data with coincident 
data collected by a conventional airborne laser altimeter. We find that freeboard estimates agree to 
within 0.03 m in the areas where sea-ice floes were interspersed with wide leads, and to within 0.07 m 
elsewhere. MABEL data may also be used to infer sea-ice thickness, and when compared with coincident 
but independent ice thickness estimates, MABEL ice thicknesses agreed to within 0.65 m or better. 
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1. INTRODUCTION 

Sea-ice mass balance is a leading indicator of the state of the 
global climate system. Understanding the nature and vari- 
ability of Arctic and Southern Ocean sea ice is critical to 
improve predictions of future climate state (NRC, 2012). 
Altimeters on both airborne and satellite platforms provide 
measurements of sea-ice elevation, from which sea-ice 
freeboard and thickness may be inferred (e.g. Laxon and 
others, 2003; Kwok and Cunningham, 2008; Farrell and 
others, 2009). By incorporating knowledge of ice concen- 
tration, sea-ice volume may be estimated. Thus the ice pack 
may be monitored remotely, and an assessment of inter- 
annual variability and decadal trends in ice thickness and 
volume may be established. 

Observations from the Ice, Cloud and land Elevation 
Satellite (ICESat), Envisat and CryoSat-2 missions indicate a 
decline in the thickness and volume of the Arctic ice pack 
over the past 10 years (e.g. Giles and others, 2008; Kwok 
and others, 2009; Kurtz and others, 201 1; Laxon and others, 
2013). The sea-ice pack has changed rapidly in its 
composition, transitioning from a predominantly thick, 
multi-year ice pack to a thinner, seasonal pack (Comiso, 
2012), and the greatest losses have been observed in the 
oldest, multi-year ice areas (Farrell and others, 2009). 
NASA's Operation IceBridge (OIB) mission (Koenig and 
others, 2010), which commenced in March 2009, includes 
yearly (every March/April) airborne campaigns over critical 
regions of the Arctic sea-ice pack, and extends the moni- 
toring provided by satellites. IceBridge provides altimetric 
measurements of the ice pack at higher resolution (~tens of 
meters) than the satellite altimeters (^hundreds to thousands 


of meters). Recent assessment of the IceBridge dataset shows 
that, following the precipitous drop in multi-year ice in 
2007/08 (Giles and others, 2008; Haas and others, 2008; 
Kwok and others, 2009), the central Arctic remains domin- 
ated by multi-year ice just over 3 m thick, while ice 
thickness in the Beaufort and Chukchi Seas appears to 
continue to decline (Richter-Menge and Farrell, 2013). 

With the planned launch of ICESat-2 in 2017 we may 
expect continuity of the sea-ice freeboard, and hence 
thickness, time series through the end of this decade. 
Accurate measurement of monthly, basin-scale sea-ice 
freeboard is a key mission requirement for ICESat-2 (http:// 
icesat.gsfc.nasa.gov/icesat2/). ICESat-2 data, when com- 
bined with observations from previous and ongoing 
missions (e.g. European Remote-sensing Satellites 1 and 2 
(ERS-1/2), Envisat, ICESat, CryoSat-2, IceBridge), will allow 
us to better understand the evolution of the ice pack over a 
multi-decadal period. ICESat-2 will carry the Advanced 
Topographic Laser Altimeter System (ATLAS), which uses a 
different detection strategy to the Geoscience Laser Altim- 
eter System (CLAS) on ICESat. ATLAS is a high-repetition- 
rate (10 kHz), multi-beam, 532 nm wavelength laser that 
utilizes photon-counting (PC) techniques for surface detec- 
tion (Abdalati and others, 2010). Specifically, ATLAS will 
have six laser beams, each with a ~1 0 m diameter footprint 
and ~0.07 cm along-track sampling, and will measure the 
range of individual photons returned from the surface. 
ATLAS beams will be arranged in three sets of beam pairs, 
with spacing between sets of pairs ~3000m apart, while 
each beam within a given pair will be ~90m apart (Brunt 
and others, 2014). This configuration is primarily designed 
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Fig. 1. Map of airborne survey over sea ice in the Greenland Sea conducted on 10 April 2012. Sea-ice conditions are illustrated with a 
contemporaneous Aqua MODIS image, acquired ~6.3 hours before the MABEL survey, and indicate landfast ice along the Greenland coast 
and large icefloes in the marginal ice zone (MIZ). Blue line shows the trajectory of the NASA P3-B, which conducted the survey, transiting from 
point (B) to (A) to (B) in an out-and-back loop, at the times indicated. Pink line illustrates the trajectory of the NASA ER-2 survey, which transited 
from point (A) to (B) in 35 min, at times indicated. Red stars indicate four case-study regions with varying sea-ice conditions. Inset shows Arctic 
region highlighting the Greenland Sea (GS), Eram Strait (ES), Beaufort Sea (BS) and Chukchi Sea (CS) with an overview of the flight surveys. The 
P3-B flight was conducted from Kangerlussuaq (SEj), Greenland, while the ER-2 flight was conducted from Keflavfk (KEE), Iceland. 


to determine surface slope in areas of high topographic relief 
on ice sheets, but will also be useful for accurate sea-ice 
freeboard determination by increasing the number of 
available elevation measurements across leads. 

ICESat-2 offers some considerable improvements com- 
pared to its predecessor ICESat, in particular the dense 
along-track sampling of the surface, which will allow 
flexibility in the post-processing algorithmic approaches 
taken to optimize the signal-to-noise ratio and obtain highly 
detailed measurements of the sea-ice pack (Farrell and 
others, 201 1). With digital PC, uncertainty arises due to not 
knowing from where within the 10 m footprint the photon 
was returned. Exploiting the dense surface sampling, along- 
track averaging is therefore employed to improve measure- 
ment accuracy (Abdalati and others, 2010). 

Here we demonstrate the utility of digital, micro-pulse PC 
laser altimetry over the polar oceans for the derivation of 
sea-ice freeboard. The Multiple Altimeter Beam Experimen- 
tal Lidar (MABEL) is an airborne simulator for ATLAS (McGill 
and others, 2013). Although the specifications for MABEL 
differ slightly from those for ICESat-2 (Kwok and others, 
2014), MABEL was designed to test the instrument theory 
and measurement strategy of ICESat-2. MABEL data are 
being used to develop geophysical algorithms prior to 
launch (Brunt and others, 2014). In this study, we assess the 
quality of the first PC altimetry dataset gathered over the 
Arctic ice pack, and its usefulness for deriving sea-ice 
freeboard. We utilize coincident, conventional analogue 


laser altimetry to determine the accuracy of the novel 
freeboard estimates derived from MABEL. Utilizing novel 
snow depth estimates from IceBridge we also assess sea-ice 
thickness estimates derived from MABEL freeboard data. 

2. DATA 

2.1 . Experiment configuration 

An airborne experiment took place on 1 0 April 201 2, during 
which the NASA ER-2 aircraft carrying MABEL, and the 
NASA P-3B aircraft carrying the IceBridge instrument suite 
(http://icebridge.gsfc.nasa.gov), surveyed sea ice in the 
Greenland Sea, south of Eram Strait (Fig. 1). Sea-ice 
conditions in the area (illustrated in a contemporaneous 
Moderate Resolution Imaging Spectroradiometer (MODIS) 
image acquired ~6.3 hours before the MABEL survey; Fig. 1 ) 
included landfast ice alongthe Greenland coast and large ice 
floes in the marginal ice zone (MIZ). The P-3B surveyed the 
region by transiting northbound (from point (B) to (A), Fig. 1 ) 
at ~460 m above the ice and returned from point (A) to (B) at 
~91 5 m, measuring a wider swath of ice on the southbound 
leg. The ER-2 aircraft overflew the P-3B at ~1 8 000 m at the 
northern end of the track, and, for operational reasons, only 
collected data on its return leg from (A) to (B). The weather 
along the survey line varied, with cloud and low-level ice fog 
north of point (A) and south of point (B), while the central part 
of the line had clear-sky conditions. These weather condi- 
tions were verified using near-coincident MODIS imagery 
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Fig. 2. Sea-ice elevation and morphology at the southern end of the survey in case-study region 4 (see Fig. 1 for location), (a) ATM elevation 
swath map with MABEL beam 6 geolocation (white line); (b) DMS imagery; (c) raw MABEL elevation profile, beam 6; (d) raw MABEL 
elevation profile, beam 5; and (e) DMS pixel brightness. Elevations are quoted with respect to the geoid, and geophysical and drift 
corrections have been applied. 


and data collected by the Cloud Physics Lidar (http://cpl.gsfc. 
nasa.gov/), a second instrument flown on the ER-2. This study 
focuses on an assessment of MABEL data quality at four 
locations (Eig. 1, red stars) between points (A) and (B), via 
comparison with data from two IceBridge instruments: the 
Airborne Topographic Mapper (ATM) laser altimeter (Krabill, 
201 3) and the Digital Mapping System (DMS) high-resolution 
camera (Dominguez, 2010). Data assessed in each study 
region (Eig. 1, red stars) span ~28km and represent a wide 
variety of sea-ice types. 

2.2. Airborne data 

The ATM airborne laser altimeter operates at a wavelength 
of 532 nm with a pulse repetition frequency of 5 kHz and a 
scan angle of ~1 5°. Depending on flight altitude during the 
survey, the P-3B mapped elevations across a swath, ranging 
from approximately 250 to 575 m. For a nominal flight 
altitude, the footprint size was ~1 m and footprint separ- 
ation at the swath center was typically ~5 m in the along- 
track direction. The ATM digitizes returned energy as a 
waveform, from which surface elevation is derived (Fig. 2a). 
Here we utilize the IceBridge ATM Level-1 B Elevation and 
Return Strength (ILATM1B) dataset (Krabill, 2013). Elevation 
measurements were interpolated onto a regular 10m grid 
using a near-neighbor routine. 


We also make use of the DMS LI B Geolocated and 
Orthorectified Images (IODMS1 B) dataset, which consists of 
high-resolution natural color imagery with a pixel resolution 
of ~0.1 m (Dominguez, 2010). The DMS images are useful 
for interpreting sea-ice features (Fig. 2b), and pixel bright- 
ness (intensity) may be used for lead discrimination (Fig. 2e). 
Here we interpolated the DMS data onto a regular 1 0 m grid. 

MABEL is a PC laser altimeter with multiple beams 
arranged in a linear array, perpendicular to the direction of 
flight, operating at two wavelengths: 1064 and 532 nm. 
MABEL beams have variable signal strengths and the system 
allows for beam geometry changes between flights if 
required. In this study we assess MABEL release 8 data from 
beams 6 and 5, which are two 532 nm wavelength beams, 
located at the center of the linear array and at ~85 m (ground 
spacing) away from the center, respectively (Brunt and 
others, 2014). We assessed beams 5 and 6 since, for this 
flight, they are most similar to the expected radiometry of the 
strong beams on ATLAS, itself a 532 nm wavelength laser 
(Brunt and others, 2014). During this experiment the pulse 
repetition frequency of the MABEL laser was 5 kHz, such that 
for a nominal aircraft speed of ~200 m s“\ a ~2 m diameter 
footprint was illuminated every ~0.04m along-track (Brunt 
and others, 2014). Raw, unfiltered MABEL data collected 
over the ice pack include signal photons reflected from the 




170 


Farrell and others: Sea-ice freeboard from photon-counting laser altimetry 


344.94' 344.96 344.98 



0.0 0.5 l.O 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

Surface elevation above geoid 


Fig- 3. Swaths of sea-ice elevation, derived from the ATM laser 
altimeter, for the northbound (bounded by dashed white line) and 
southbound legs of the IceBridge sea-ice survey on 10 April 2012. 
Local ice-drift vectors (black arrows) were estimated by manually 
tracking the geolocation of features in the sea-ice pack (e.g. red 
circles indicate a deformed sea-ice floe). 


sea-ice/ocean surface, and background photons due to solar 
backscatter in the atmosphere and detector noise. For this 
flight ~4.3 signal photons per meter along-track were 
recorded for beams 5 and 6 (Brunt and others, 2014). An 
absolute elevation bias correction was applied to the MABEL 
data, following Brunt and others (2014), to account for 
differences in the optical paths, unique to each MABEL 
beam. For beam 5 the bias correction was 1.99 m, while it 
was 1 .76 m for beam 6. Figure 2c and d show MABEL along- 
track measurements for beams 6 and 5, respectively, and 
indicate signal photons close to the surface and background 
photons dispersed in the atmospheric column. 

2.3. Geophysical corrections 

The Earth Gravitational Model 2008 (EGM2008; Pavlis and 
others, 2008) was subtracted from all altimetric elevation 
measurements (ATM and MABEL) to provide surface height 
above the geoid. The GOT 4.8 ocean tide correction was 
applied to remove the effects of local tides during the time 
frame of the experiment and includes a long-period equi- 
librium ocean tide correction to remove tidal effects due to 
the Sun (Ray, 1 999). 


2.4. Sea-ice drift correction 

Since each NASA aircraft operated at different flight speeds, 
there was a shift in temporal sampling between the two 
airborne surveys. High rates of sea-ice drift are encountered 
south of Eram Strait and this complicates assessment of data 
gathered by the two aircraft systems in this region. To 
achieve accurate intercomparisons we must therefore ac- 
count for ice velocity. We estimate and correct for the local 
sea-ice drift rate by comparing the geolocation of specific 
sea-ice features (e.g. pressure ridges and rubble fields) in the 
northbound and southbound ATM elevation swaths (Fig. 3). 
Velocity vectors were computed to correct for the temporal 
offset between the IceBridge and MABEL data. We estimate 
an average drift rate of 0.25 ms“\ which agrees to within 
0.01 m s“^ with an independent estimate of sea-ice drift on 
the date of the experiment, derived from a combined 
Advanced Scatterometer and Special Sensor Microwave/ 
Imager product available at CERSAT/lfremer (http://cersat. 
ifremer.fr/). Combining knowledge of the drift rate, and the 
time elapsed between MABEL and IceBridge acquisitions, 
we adjusted the geolocation of the MABEL data to coincide 
with the ATM and DMS data. 

3. METHODS 

We apply the MABEL surface-finding algorithm outlined in 
Brunt and others (2014), a histogram-based approach for 
distinguishing signal photons from background photons, to 
reduce the MABEL point cloud such that it only includes the 
coarse signal photons associated with sea-ice and ocean 
surface returns (Eig. 4a). We apply a one-dimensional (1-D) 
along-track maximum likelihood probability filter, operating 
at 10 m, to the MABEL data to further reduce noise and 
smooth the surface response of the lidar (Fig. 4b) to achieve 
consistency in the spatial sampling of both the MABEL and 
IceBridge datasets. 

Data from the laser altimeters assessed in this study 
provide measurements of both sea-ice and lead elevations. 
Eor ice floes, the primary reflection is from the air/snow 
interface via surface scattering. Sea-ice freeboard, the 
portion of an ice floe above local sea level, is therefore 
defined here as the ice-plus-snow freeboard. Accurate 
identification of leads between floes is required to obtain 
measurements of the sea-surface height (SSH, also known as 
'tie points'). DMS visible imagery, spatially coincident with 
both MABEL and ATM, was assessed. Normalized pixel 
brightness (0-255) at each point along the lidar beam-path 
was extracted (Fig. 2e). Using a pixel-brightness threshold 
range of 0-50 to classify leads, we flag signal photon 
reflections from lead surfaces (Fig. 4b, light-green dots). We 
identify ATM lead elevation measurements in the same way 
(Fig. 4c, orange dots). The threshold value for lead classifi- 
cation was derived empirically, based on close visual 
inspection of the imagery, and was chosen to include open 
water, nilas and thin grey ice, while excluding refrozen 
leads with snow accumulation. Here, modal lead elevation 
defines the local SSH and was linearly interpolated between 
tie points to estimate along-track SSH in each 28 km study 
region (Fig. 4b and c, red lines). Along-track SSH is 
subtracted from sea-ice elevation measurements to derive 
sea-ice freeboard. 

To demonstrate the denser along-track surface sampling 
afforded by digital PC altimetry, we indicate a set of simu- 
lated surface elevations assuming the ICESat configuration 
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Fig. 4. Sea-ice elevation measured along-track in case-study region 4. (a) MABEL elevations from beam 6, after application of surface-finding 
algorithm, (b) Same as for (a), after application of 10m along-track filter, (c) ATM elevation profile. Lighter-color dots (light-green (b) and 
orange (c)) indicate the location, and height, of lead elevation measurements, after application of lead discrimination threshold derived from 
DMS imagery (Fig. 2e). Red line (b, c) indicates sea surface height, interpolated along-track. Open red circles (a) indicate a set of simulated 
surface elevations assuming the ICESat sampling configuration. Segment shown is ~28km long. 


(Fig. 4a, open red circles). The GLAS on ICESat consisted of a 
single-beam laser altimeter with a ~50-70m footprint 
diameter, spaced every ~172m along-track (Farrell and 
others, 2009). Following the approach outlined in Farrell and 
others (2011) we sampled the MABEL surface elevations 
using the ICESat/GLAS footprint pattern, recording the mean 
elevation per footprint using a 50 m diameter search radius. 
Application of the MABEL surface-finding algorithm (Brunt 
and others, 2014) across the 28 km long segments assessed 
here resulted in ~75 000 surface elevation measurements, 
which reduced to ~3200 measurements after application of 
the 1-D 10m along-track filter. The ICESat sampling ap- 
proach would have provided only ~1 70 elevation measure- 
ments across the same distance. 

3.1. Measurement accuracy 

The single-shot accuracy of the ATM elevation measure- 
ments is ~0.1 m (Krabill and others, 2002), and geolocation 
accuracy is estimated at 1 m (Schenk and others, 1999). 
Interpolation of ATM data onto a regular grid results in an 
elevation precision of ~0.05 m over level, undeformed sea 
ice (Earrell and others, 2012). The accuracy of the MABEL 
elevation measurement is affected by geolocation-induced 
elevation uncertainty, time-tag errors within MABEL, and 
errors in the surface-finding algorithm. For sea-ice floes with 
low surface slope, MABEL elevation uncertainty is estimated 
at ~0.15m (Brunt and others, 2014). Brunt and others 
(2014) note that while MABEL elevation uncertainty is larger 
than that of other airborne lidars, it is within the instrument 
design specifications for ATLAS algorithm development. 
Assessing the standard deviation of elevation measurements 
over nilas and new, level sea ice in leads, we also determine 


a MABEL elevation precision of 0.1 5 m, consistent with the 
results of Brunt and others (2014). This further reduces to 
0.04 m, after application of the 10m along-track filter. 

4. RESULTS 

We assessed sea-ice elevation measurements acquired by 
both the MABEL and ATM instruments at four case-study 
locations along the airborne survey line (Fig. 1, red stars). 
These four ~28km-long segments were associated with a 
wide range of sea-ice conditions, outlined in Table 1 . After 
filtering, the MABEL elevation data reveal fine details of the 
surface roughness and pressure ridging of Arctic sea ice 
undergoing export through the Fram Strait, while leads and 
nilas are also well defined (e.g. Fig. 4b). We find strong 
agreement between the sea-ice morphological features 
identified in the MABEL data (Fig. 4b) and those revealed 
in the ATM surface elevation measurements (Fig. 4c) and in 
the DMS imagery (Fig. 2b). 

We derived sea-ice freeboard distributions for case-study 
regions 1, 3 and 4, comparing freeboard estimates derived 
from MABEL beams 5 and 6 with freeboard derived from the 
ATM data. Study area 2 did not contain any leads except for 
sparse fractures. Thus neither freeboard nor thickness could 
be determined in case-study region 2, and statistics de- 
scribing the elevation distribution are provided instead. 
Freeboard distributions for regions 1, 3 and 4 are shown in 
Figure 5. Mean and modal freeboard statistics are presented 
in Table 1 . Overall we found good agreement in the form of 
the derived freeboard distributions for both MABEL and ATM 
data in all study regions. The bimodal distributions represent 
the SSH (lower mode) and the sea-ice freeboard (upper 
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Fig. 5. Sea-ice freeboard distributions for case-study regions 1 (a), 3 (b) and 4 (c). Results are shown for MABEL beam 6 (green), beam 5 
(blue) and IceBridge ATM (magenta). Bin size is 0.03 m. 


mode), while the elongated tails indicate the presence of 
thick, deformed ice floes in the survey area. Mean MABEL- 
derived freeboard agreed with mean ATM-derived freeboard 
to ±0.07 m across all study regions. The best agreement 
(±0.03 m) was found in study region 3, which was composed 
of discrete floes with pressure ridges and rubble, interspersed 
with wide leads. Modal freeboard values were typically 
lower than mean freeboard values and varied by 0-0.15 m 
across instruments and study regions. In the region where no 
leads were present (case study 2), we find that the application 
of geophysical range corrections (Section 2.3) and the 
application of the MABEL bias correction for beams 5 and 
6 (Section 2.2) provided mean surface elevations consistent 
to within 0.09 m or better. 

Combining the derived MABEL and ATM freeboards with 
novel snow depth estimates available from IceBridge, we 
estimate the thickness of sea ice in case-study regions 1, 3 


and 4. Thickness may be inferred from altimetric measure- 
ments by assuming hydrostatic equilibrium, estimates for the 
density of snow, ice and sea water, and combining ice 
freeboard with knowledge of snow depth. EHere we calculate 
mean sea-ice thickness following Giles and others (2007, 
eqn 2). Separately, data from the full IceBridge instrument 
suite are processed and released within 6 months of each 
IceBridge field campaign. These data are used to derive sea- 
ice products including freeboard, thickness and snow depth 
(Kurtz and others, 2013a) that are publicly available via the 
US National Snow and Ice Data Center (NSIDC). Elere we 
assess the quality of sea-ice thicknesses derived from the 
retrieved MABEL and ATM freeboards via a comparison with 
the independently derived IceBridge thickness product 
(IDCSI2; Kurtz and others, 2012) for 10 April 2012. In order 
to make an appropriate comparison we follow the values 
used by Kurtz and others (201 3a, and references therein), for 


Table 1. Mean sea-ice freeboard and thickness estimates derived from MABEL and ATM elevation measurements for each study region. 
Modal freeboard results are based on the distributions provided in Figure 5 using a bin size of 0.03 m. A description of sea-ice morphology 
and atmospheric conditions in each study region is provided. Snow and sea-ice thickness provided in the IceBridge sea-ice product are 
included for comparison with estimates from MABEL and ATM 




1 

Case-study 

2 

area 

3 

4 

Sea-ice conditions 


Ice floes with pressure 

Consolidated floes with 

Discrete floes with 

Discrete floes with 



ridges and narrow 

rubble, fast ice, sparse 

pressure ridges and 

pressure ridges and 



refrozen leads 

narrow fractures, no leads 

rubble, interspersed 

rubble, interspersed 





with wide leads 

with wide leads 

Atmospheric conditions 


Clear 

Clear 

Clear 

Clear 

Mean sea-ice freeboard* 

ATM 

0.55 

0.89 

0.72 

0.62 

(m) 

MABEL beam 6 

0.57 

0.95 

0.75 

0.58 


MABEL beam 5 

0.6 

0.98 

0.73 

0.55 


OIB product 

0.53 

0.85 

0.73 

0.57 

Modal sea-ice freeboard* 

ATM 

0.475 ±0.01 5 

0.825 ±0.015 

0.750 ±0.01 5 

0.475 ±0.01 5 

(m) 

MABEL beam 6 

0.625 ±0.01 5 

0.925 ±0.01 5 

0.700 ±0.01 5 

0.625 ±0.01 5 


MABEL beam 5 

0.575 ±0.01 5 

0.875 ±0.01 5 

0.725 ±0.01 5 

0.625 ±0.01 5 

Mean snow thickness (m) 

OIB product 

0.24 

0.23 

0.28 

0.24 

Mean sea-ice thickness 

ATM 

3.60 (0.60) 


4.94 (0.60) 

4.29 (0.60) 

(uncertainty) (m) 

MABEL beam 6 

3.79 (0.53) 


5.22 (0.53) 

3.92 (0.53) 


MABEL beam 5 

4.07 (0.53) 


5.04 (0.53) 

3.64 (0.53) 


OIB product 

3.97 (0.61) 

6.09 (0.76) 

5.29 (0.71) 

4.15 (0.63) 


*Neither sea-ice freeboard nor thickness could be estimated in study region 2 due to a lack of leads. Surface elevation statistics are provided instead. 
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the density of snow, ice and sea water. Snow thickness 
(Table 1) is taken directly from the IceBridge product (Kurtz 
and others, 2013a). Mean sea-ice thickness for ATM and 
MABEL beams 5 and 6 is provided in Table 1 and can be 
compared to the independently derived mean ice thickness 
estimate provided in the IceBridge product for the same study 
region. The sea-ice thickness estimates from MABEL and 
ATM are on average within 0.37 m of each other and these 
are in turn consistent with the IceBridge product. Average 
sea-ice thickness in the Greenland Sea study area is 4.25 m, 
and independent thickness estimates differ by 0.07-0.65 m. 

Sea-ice thickness uncertainty estimates are also provided 
in Table 1 . Following Giles and others (2007, eqn 4) an error 
propagation analysis is utilized to estimate thickness 
uncertainty assuming the following terms: The uncertainty 
associated with the IceBridge snow thickness is 0.057m as 
provided in the product and described by Kurtz and others 
(2013a). The uncertainty associated with an individual ATM 
freeboard estimate is equated to the measurement precision 
of ~0.05 m for interpolated ATM data over level, unde- 
formed sea ice. Likewise the uncertainty associated with an 
individual MABEL freeboard estimate is equated to the 
measurement precision of ~0.04m for MABEL data over 
refrozen leads after the application of the 10 m along-track 
filter (Section 3.1 ). The differences between independent ice 
thickness estimates fall within the uncertainty for each 
instrument/product, which ranges from 0.53 to 0.76 m. 


5. SUMMARY AND CONCLUSIONS 

We have demonstrated that MABEL digital, micro-pulse PC 
laser altimetry data are effective for reliably determining the 
freeboard and thickness of Arctic sea-ice floes, when 
combined with auxiliary information on snow depth and 
snow/ice density. The surface-finding algorithm efficiently 
distinguished coarse signal from background photons to 
allow assessment of elevation measurements associated with 
both ice floes and leads. The surface finder tracked the 
surface of both rough sea-ice floes and smoother lead 
surfaces, providing unprecedented details of the surface 
morphology of Arctic sea ice. We have verified the accuracy 
of the MABEL freeboard estimates through comparison with 
independent, coincident measurements collected by the 
ATM laser altimeter on board IceBridge. The results show 
that the surface elevation and freeboard distributions derived 
from the two independent airborne lidar sensors (MABEL and 
ATM) are consistent. The ATM- and MABEL-derived sea-ice 
freeboard estimates agreed to within 0.07 m or better across 
all study areas, and consistent freeboard estimates (0.03 m or 
better) were derived in study region 3 where deformed ice 
floes were separated by wide leads. We have also demon- 
strated that MABEL data may be used to successfully infer sea- 
ice thickness. MABEL ice thickness estimates when com- 
pared with coincident data from the ATM and the IceBridge 
thickness product agreed to within 0.65 m or better. 

Although sea-ice freeboard is a differential measurement, 
knowledge of absolute range is required to define the 
precision and overall quality of PC elevation measurements 
of the ice pack. Measurements across undeformed areas of 
sea ice (e.g. nilas and new ice in leads, level ice pans or level 
snowfields on fast ice) could be used to calibrate the laser 
signal, calculate measurement precision, assess the impact of 
along-track sampling and filtering and investigate potential 
geophysical range biases due to volume scattering at the 


surface. In this study, the assessment of MABEL digital PC 
data quality and accuracy was greatly enhanced by the 'out 
and back' flight-loop conducted by the IceBridge aircraft. 
This survey configuration constrained the local ice-drift 
vectors and allowed rectification of the two airborne datasets 
to a common reference time frame to mitigate the temporal 
offset between MABEL and IceBridge acquisitions. The 
airborne experiment conducted on 10 April 2012 serves as 
an excellent template for future MABEL missions aimed at 
assessment of photon-counting laser altimeters over sea ice 
and provides a demonstration of techniques that could be 
used to validate ATLAS measurements collected by ICESat-2 
over the global sea-ice packs after launch. To assess MABEL/ 
ATLAS digital PC data, a second altimeter with sufficient 
accuracy and across-track spatial sampling (i.e. scanning/ 
swath mapping) is desirable so as to capture the multi-beam 
MABEL and ATLAS beam configurations. In the case of 
ATLAS, the central beam pair will be separated by ~90m, 
while two outer beam pairs will each be ~ 3000 m from the 
central beam pair (fig. 1 in Brunt and others, 2014). 

Further work will include an assessment of additional 
MABEL data gathered over sea ice during the April 2012 
flight campaign and data from future sampling campaigns 
aimed at obtaining data under different seasonal and 
atmospheric conditions. The goal for future MABEL sea- 
ice deployments is to collect measurements that will 
support the development of operational algorithms for the 
retrieval of sea-ice geophysical parameters from ATLAS on 
ICESat-2 (Kwok and others, 2014). Following the analysis 
presented in this study, we will investigate an automated 
approach for deriving sea-ice freeboard estimates from 
MABEL data that are specifically designed for use in the 
MIZ. This would take advantage of measurements across 
numerous lead surfaces that provide plentiful SSH tie points. 
We will explore the feasibility of producing a set of 'quick 
look' sea-ice freeboard estimates from ICESat-2 measure- 
ments in the MIZ, similar to the quick-look sea-ice products 
provided as part of the IceBridge mission (Kurtz and others, 
2013b). This could provide information on sea-ice free- 
board and thickness in an area of primary interest to the 
operational sea-ice community. 
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